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ABSTRACT 

Growth rates for gravitational instabilities in a thick disk of gas and stars are 
determined for a turbulent gas that dissipates on the local crossing time. The 
scale heights are derived from vertical equilibrium. The accuracy of the usual 
thickness correction, (1 + kH)~^, is better than 6% in the growth rate when 
compared to exact integrations for the gravitational acceleration in the disk. 
Gas dissipation extends the instability to small scales, removing the minimum 
Jeans length. This makes infinitesimally thin disks unstable for all Toomre-Q 
values, and reasonably thick disks stable at high Q primarily because of thickness 
effects. The conventional gas+star threshold, Qtot, increases from ~ 1 without 
dissipation to 2 or 3 when dissipation has a rate equal to the crossing rate over a 
perturbation scale. Observations of Qtot ~ 2 — 3 and the presence of supersonic 
turbulence suggest that disks are unstable over a wide range of scales. Such 
instabilities drive spiral structure if there is shear and clumpy structure if shear 
is weak; they may dominate the generation of turbulence. Feedback regulation 
of Qtot is complex because the stellar component does not cool; the range of 
spiral strengths from multiple arm to flocculent galaxies suggests that feedback 
is weak. Gravitational instabilities may have a connection to star formation even 
when the star formation rate scales directly with the molecular mass because the 
instabilities return dispersed gas to molecular clouds and complete the cycle of 
cloud formation and destruction. The mass flow to dense clouds by instabilities 
can be 10 times larger than the star formation rate. 



Subject headings: Instabilities — Stars: formation — ISM: general — Galaxies: 
structure 
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Introduction 



The stabilit y of a galactic dis k with gas and stars was first considered i n the context of 
stellar waves by iLin &: Shul (Il966[ ) and in the context of cloud formation by |Jog &: Solomon 
(11984). Th e impo rtance of cold gas in the general stability of a combined disk was recognized 
by iToomrd (Il964[). who also considered the effects of thickness and particle dynamics for the 
stars, as did iKalnajsl (119651 ) . Two-compo nent d i sks o f fini te thickness were consid ered in 
more detail for the spiral wave pr o blem by Katd (119721) and lBertin fc Romed ( 1l988l ) and in 
the general case by iRomed ( 1l992l ). iRafikovl (120011 ) investigated multiple components. 



With these two-component analyses came stability parameters that de termined in ab- 
solute terms when a disk was stable. For one stellar component, this is the iToomrel (Il964[ ) 
value for the minimum velocity dispersion in the radial direction, written now as a dimen- 



sionless parameter Q. iGoldreich &: Lynden-Belll ( 11965! ) evaluated a maxin aum density for 



the st a bility of a one-com pon ent disk. F ollowing the two-fluid analysis by jjog &: Solomon 
( ll984j ) , lElmegreenI ( 11995! ) and iJogl (Il996[ ) deter nained a comb i ned Q for a two-fluid disk. A 
sir npler form for the combine d Q was w r itten by IWang fc Silk! (!l994j ) , and recently improved 
by !Romeo fc WiegertI (120111 ). !Rafikov! (12001! ) discussed disk stability using the dispersion 
relation at zero frequency, which is a function of wavenumber, and gave an expre ssion that 
was compared to observations for spiral galaxies and dwarfs by lLeroy et al.! (120081 ). Leroy et 
al. evaluated the Rafikov expression assuming a range of likely wavenumbers, and showed 
that most radii in all of the galaxies they studied were marginally stable. This common 
feature suggests that galaxies regulate their stability by varying the velocity dispersions of 
gas and stars through star formation feedback and stellar wave scattering. 

An important component of disk stability that has been overlooked in these studies 
is the dissipation of energy in the interstellar medium. Usually the gas is assumed to be 
isothermal or polytropic, which has a certain level of dissipation by assumption, but the 
nature of this dissipation is unrealistic as it implies an ever-present source of energy to 
accompany any decompression. On average, lower-density gas is warmer than higher-density 
gas, but this relationship is established in the real ISM by specific energy sources and sinks 
that do not automatically heat up the gas if it is driven to a lower density by gravity. A 
time-dependent energy equation should be used instead of an adiabatic pressure- density 
relation. In equilibrium, this energy equation states that heating balances cooling, but out 
of equilibrium, when a spiral arm or cloud begins to form in an instability, there is usually a 
net cooling from dissipation in turbulence. Only after the spiral or cloud make stars will the 
energy loss be returned, but this is long after the instability has grown. Recent computer 
simulations suggest that the primary role of star formation feedback in a galaxy model is to 
break apart the dense clouds that form; turbulence on larger scales is driven by gravitational 
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instabilities (Bournaud et al. 2010, Hopkins et al. 2011, see also Elmegreen et al. 2003). 
This means that turbulent dissipation accompanies the growth of gravitational instabilities. 

The purpose of this paper is to introduce the energy equation to the interstellar com- 
ponent of the stability problem for a gas and star disk. We also consider thickness effects 
in some detail, checking the approximation often used for stability analyses with an exact 
integration for the gravitational acceleration in a three dimensional disk. Then we discuss 
corrections to the commonly used stability thresholds that arise from gaseous energy dissi- 
pation. The main point is that energy dissipation changes the dispersion relation for disk 
instabilities in a fundamental way, removing the Jeans length from the problem in some 
cases, and giving residual instabilities at high wavenumbers even when the Toomre Q pa- 
rameters are large. Disk thickness is a more effective stabilizer than interstellar turbulent 
pressure. 

The results have implications for the interpretation of observations, primarily in sug- 
gesting a better stability threshold when dissipation is important, and they have implications 
for simulations of galaxy disks. In simulations, an equation of state for the gas has to be 
considered to connect density changes with pressure changes. If this equation of state is 
too simple, particularly if it does not reproduce the expected rapid dissipation of turbulent 
energy in perturbations, then the simulation is probably more stable than a real galaxy - at 
least on small scales. 

Another application is to regions of galaxies where the usual stability parameters sug- 
gest a high degree of stabihty, such as the far-outer regions of d i sks, or in dwarf ga laxies 



(IHunter fc Plummerlll996l : iMeurer et al.lll996l : Ivan Zee et al.lll997l : iHunter et al.lll998l ). The 



present results suggest that these regions can still be unstable at a low level as perturbations 
in density lead to a simultaneous increase in self-gravity and a decrease in gaseous velocity 
dispersion because of dissipation. Such a region will collapse if dissipation operates on about 
a crossing time, because the pressure does not increase fast enough to overcome gravity. 

In what follows, the star-f-gas dispersion relation is derived in Section [2l a simple example 
including turbulent dissipation but without stars is discussed in Section El and the star-|-gas 
relation is evaluated in Section HJ Gas and star growth rates are considered separately 
in Section 14.11 and together in Section 14.21 The ratio of perturbed stellar and gaseous 
column densities is in Section [^^31 Section discusses dissipation corrections to the popular 
gas+star Q thresholds that are in the literature. Section |5] examines thresholds for azimuthal 
instabilities. A discussion of several implications of this work is in Section El including 
the greater role played by instabihties when dissipation is included (Sect. 16. ip . feedback 
regulation of disk stability (Sect. 16. 2p . and star formation (Sect. 16. 3p . A summary is in 
Section [71 
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The scale height of the disk enters the thickness correction factor in the gravitational 
forcing term. Scale heights are evaluated in a self-consistent way in Appendix A and used for 
the solutions of the dispersion relation in Section HJ The accuracy of the thickness correction 
factor is determined in Appendix B. 



2. Dispersion Relation 



The equations f or perturbed r adial and azimuthal motions of gas in an axisymmetric 
galaxy disk are (e.g.. lRafikovll200l[ ) 



2Qvp-o = — -— 



1 dP 



dr pgfi dr ' 



dv„e 



2Bv„, = 0. 



(2) 



dt ^ 

Here, fg^r and Vg^o are the perturbed radial and azimuthal gas velocities, is the perturbed 
gravitational potential, P is the perturbed pressure and is the unperturbed gas density. 
The angular rotation rate is Q, and B is the Oort constant of rotation. 

Perturbations of the form e*^^^"^'*) are typically introduced. The oscillation frequency 
is denoted by to distinguish it from a growth rate u introduced later, i.e., u = —iu^. The 
equations of motion are combined with the definition of the epicyclic frequency, k'^ = —4QB, 
to give 



— UJ? 



(3) 



The stars were originally treated by iToomrd (119641 ). iKalnajd (119651 ). and iLin fc Shu 



(11966! ) using the coUisionless Boltzmann equation for the phase-s pace distribut ion function. 
After analogous reduction, the stellar equation becomes (see also lRafikovll200ll ) 



H? — bj} 



where 



sin(7rsi) 



exp (— ^^[l + cos^]) sin(si^) sin^rf^ 



(4) 



(5) 



(IBinney &: Tremaind l2008l ). The notation here is that s\ is the dimensionless oscillation 
frequency, uj\/ to distinguish it from a dimensionless growth rate, s = —isi, introduced 
below. The radial velocity dispersion of stars in a presumed Schwarzschild distribution 
function is as', we define dimensionless wavenumber q = kag/n- 
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The left-hand side of Figure 1 shows J^{si, q^) as a function of dimensionless squared 
wavenumber for various dimensionless oscillation frequencies Sj. When si is small, 



Hsuq')^ , . (6) 



1 - + 

In the figure, the solid lines of various colors are the exact v alues of and the dashed lines 



with the corresponding colors are the approximate values. iToomrd ( 119641 ) recognized this 
similarity for Si = 0. The deviation between the two functions is too small to be seen in most 
of the figure. A blow-up insert shows how similar the solutions are for Si = 0.1 (green curves). 
Si = (blue curves), and s\ = 0.5 (magenta curves). The similarity between the exact and 
approximate curves for small Sj means that the instability condition for a collisionless stellar 
fluid, which is derived from the dispersion relation in the limit when Si = 0, is nearly the 
same as the instability condition for a fluid having the same velocity dispersion and mass 
column density. This is why the condition Qs = hcTs/ (S.SGGSs^q) < 1 for stars has the 
constant 3.36 in it, while the condition Qg = kg^/ (vrGSg^o) < 1 for gas has the constant tt. 
This constant, 3.35828, is the maximum value of 27rJ-'(0, g^)/g over variations in q. 

Equations (E]) and (jl]) contain the radial velocity, gravitational potential, and pressure. 
The latter two are related to the mass column density. This column density is related to 
velocity through the continuity equation, which is, for perturbed variables, 

WiSg = /cfg^rSg^o- (7) 

WiSs = A;t;s,rSs,o- (8) 

Column density is related to pressure through a time-dependent energy equation. The 
equation for perturbed energy, when multiplied by 7 — 1 for adiabatic index 7, becomes an 
equation for perturbed pressure changes: 

The perturbed rate of increase of energy per unit volume is F and the perturbed rate of 
decrease from radiative cooling and other effects is A. We consider that star formation is 
mostly triggered by the disk instabilities that we want to study, so that before star formation 
begins, there is very little perturbation in the heating rate. Perturbed heating might come 
from background supernovae, for example, and depend on the gas density so that the rate 
changes when the density changes. However, the heatin g rate from supe rnovae depends 



very weakly on density (as p""'^, as may be derived from lCioffi et al.l 119881 ). so we set the 
perturbed F equal to zero. This means that the background turbulent heating rate per unit 
volume does not increase much as the density increases during the growth of a perturbation. 
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Cooling comes mostly from shocks in an interstellar medium that is supersonically tur- 
bulent. The cooling rate should be larger where the kinetic energy density is larger, which is 
where the pressure is also larger. To the extent that we can approximate the momentum flux 
for a supersonically turbulent fluid with the pressure term, P, we can also approximate the 
dissipation time for the corresponding energy density as some number, 6~^, of crossing times 
in the perturbed region. Thi s crossing time dep endence follows from numerical simulations 



m 



Mac Low et al.l (Il998l ) and lStone et al.l (Il998l ). Because the crossing rate is a^k, we write 



the dissipation rate as 

(7-l)(-A) = -5agfcP. (10) 

It is significant that the perturbed cooling rate peaks at the same position as the perturbed 
pressure, rather than at the position of peak pressure gradient. Then the time dependence 
of the pressure in equation IQ is in phase with the pressure itself. 



^ \ Pgfl J \dt 

This phase coherence means that the pressure has a component to the time dependence 
that is exponential. There are no purely oscillatory solutions as there are for polytropic or 
isothermal gasses. Thus we need solutions of the form e*^'^"'"'^* for real u. We also expect that 
for self-gravitating media there will be growing (unstable) solutions, in which case a; > 0. 
Note that the derivation of J-" also requires the real p art of the growth term to be positive 



so that an integral converges at negative infinite time ( iBinney fc Tremaind 120081 . Appendix 



6A). We identify u as the growth rate of the gas+star gravitational instability. 

The reduction factor includes solutions where the time dependence is exponential. 
Substituting is for Si in equation (j5]), we obtain for growing perturbations, 

^growth(s, g^) = -T-rj — T / exp (-q^[l + cos^]) sinh(s^) sin^d^ (12) 
smh(7rs) Jq ^ 

where s is the dimensionless growth rate, s = u/k. The right-hand side of Figure 1 shows 
J-" for growing solutions. An excellent approximation is equation (j6]), which is now written 



1 + 

J^growth(g,g') - i^g2'^g2 - (13) 



This approximation, shown by the dashed lines, is closer to the exact growth solution than 
equation ([6]) was to the oscillatory solution, equation (jS]). 

The dissipation term in equation iQ has another important effect because it removes 
the minimum unstable length, or Jeans' length, from the problem for an infinitesimally thin 
disk. When 6^0, arbitrarily small perturbation lengths can dissipate turbulent energy 
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and drive an instability. It also removes a stability condition: gas perturbations that always 
cool are always unstable. Then there is no minimum gas column density for instability, 
nor is there anything analogous to a Toomre Q condition for an infinitesimally thin disk. 
Note that zero-order heating and cooling are assumed to be in equilibrium. The turbulent 
dissipation of zero-order pressure is balanced by energy input from background sources, such 
as supernovae. But to the extent that perturbed pressures driven by self-gravity always 
dissipate internal turbulent energy faster than they gain turbulent energy from background 
sources, there is an instability at all length scales. In practice, that means all length scales 
where the gas is supersonically turbulent. We return to this point in Section H] where the 
effects of disk thickness lead to a stability condition again. 

At this point we introduce perturbation variables that clearly indicate their spatial phase 
dependence. We let perturbed pressure, density, column density and gravitational potential 
vary as cos kr, and we let both velocities vary as sin kr. All perturbed quantities varying in 
time as e'^*. It follows that equations ([3]) and (jl]) become 



ku 



P/P,,o) , (14) 



Vs,r = ^^^^^growth(w//t, [kajnf). (15) 
The continuity equations are 

wSg = —kVg^r^gfi, (16) 

= -kVs,r^s,0- (17) 

The energy (pressure) equation simplifies to 

7^0 Sg f u 



P 



Pressure goes to zero as u goes to zero if 5 7^ 0. This is why pressure does not enter the 
stability condition and why there is no minimum length in the case of turbulent dissipation 
like that given by equation ffTOj) . 

The final link is between column density and gravitational potential. This comes from 
Poisson's equation for the gravitational potential and has no time dependence. The equation 
is three-dimensional, so the perturbed potential depends on the distribution of perturbed 
density perpendicular to the plane. This leads to correctio n factors for th e mass column 



densities if only the two-dimensional forcing is considered. iToomrd (119641 ) used a correc 



tion factor o f (1 — e ^^)/{kH) for a one-component uniform disk between z = ±if, while 



VandervoortI (jl970l ) used (1 + kH) ^ from a more detailed analysis. These two factors are 
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similar, as is the J factor in Shu (1968, see Vandervoort 19 70). The use o f analogous cor- 



rection factors in the two- component case was discussed by iRomed (Il992[ ). He concluded 
that the correction factors cannot generally be applied independently to each component, 
but if the vertical distributions are somewhat uncoupled, then that approximation should 
be all right. In this case the result is a perturbed gravitational potential that depends on 
the perturbed gaseous and stellar column densities with separate and independent correction 
factors, 

ItxG / So- Ss \ , ^ 

+ • (19) 



k \l + kHg l + kH,y 
We use this approximation here. In Appendix |X1 the two scale heights H are evaluated from 
the equations of vertical equilibrium. In Appendix |Bl we check the thickness approximation 
in equation ( fT9|) by direct integration over the combined density distribution to determine 
the gravitational acceleration from a perturbation. The result suggests that the (1 + kH)~^ 
correction factors for column density are too low by ~ 12% at kH ~ 1 — 10, but this is the 
largest error introduced by the thickness approximation for a wide range of Sg o/Ss,o and 
other parameters. 

With equations f|T5l) . f|T71) . and f|T9l) we can eliminate (p and Vg^r to get an equation for 
perturbed Eg/Eg: 

(1 + kH^) _ 



E^{l + kH,) l-Fs' 

where 



(20) 



_ 27rG'A;Es,oJgrowth(-s,g^) 
^~ iu^ + K^)il + kH,) ^ ' 

and s = u/ n, q = kcs/ k, as defined above. 

Now equations ([HD, ffTB]) . ( ITS]) . (IT^ . and ( 12U]) combine to give the dispersion relation, 

1 = + y; (22) 

where 

Y = 27iGkE^^o 

and the square of the unperturbed gas velocity dispersion is = 7-Po/Pg,o- With dimen- 
sionless variables, this dispersion relation is 

1 _ ^ 2gJgrowth(g, 9^) 1^24) 



Qg (1 + s2 + q^S^s/[s + 6qS]) (1 + gi^g) Qs (1 + s^) (l + qH^) 
where s = u/ k, q = kcr^/ n, S = o-g/os, Hg = Hgn/cTs, = H^n/a^, Qp = Ka g /fTrC S p. n) and 



Qs = KCTs/ (nGT^sfi). This expression is analogous to that in I Jog fc SolomonI (Il984[ ). iRomeo 
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(1992 


) and 


Rafikov 


(2001) 



our use of growth solutions for the reduction factor T . 



In iRafikovl (l200ll ). the instabihty condition for a gas+star fluid was that the right-hand 
side of equation (|24|) . written with 5 = and iJg = i/g = in his case, had to exceed 1 in 
the limit of s = for some value of wavenumber q. Stability requires the right hand side to 
be less than 1 at s = for all wavenumbers. If 5 = 0, we have these conditions too, because 
then the gas term varies as qSj (1 + g^S*^) which has a maximum value of 0.5 at qS = 1, and 
then sufficiently small Qg makes the right hand side larger than 1. When S ^ 0, however, 
the gas term varies directly with q when s = and there is no maximum value. The right 
hand side can exceed 1 for s = even for very large Qg, which is normally stable, if q is 
large enough. For an infinitesimally thin disk, we need q > Qg/2S or k > fi;^/(27rG'Sg o) for 
instability in the pure gas case (Sg q = 0) if there is perturbed gaseous energy dissipation of 
any magnitude, i.e., 6 > 0. This is the same as the instability condition for a infinitesimally 
thin pressure-less fiuid, as derived by Toomre (1964, eqn. 21). 



3. Jeans Relation with Dissipation 

The role of turbulent dissipation in the standard Jeans analysis can be seen from simpler 
equations involving only the gas. Equation (1^2]) with the first term on the right is for gas only, 
but that also contains the effects of rotation and finite thickness. Solutions to this equation 
will be shown in the next section. An even simpler form is to write it for an infinitesimally 
thin disk without rotation in the form 

tu^ = 2-KGT.gk - . (25) 

For the present discussion, we make this equation dimensionless by defining the unit of 
frequency-squared to be 27rG'Sfcj and the unit of wavenumber- squared to be fcj = 2TiGT,g/ a^; 
kj is the conventional Jeans wavenumber. The two-dimensional Jeans equation with this 
normalization is 

u;^ = k- (26) 



u + 6k 

This is a cubic equation. 



uj^ + Lo^k5 + co{k^-k)-k'^S = 0, (27) 



that can be solved analytically, although the steps are complicated with terms in 6 up to 
order 6. 
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Figure |2] shows numerical evaluations of this analytical solution as solid lines for 6 = 0, 
0.5 and 1. The usual dispersion relation with a wavenumber limit for growth at k/kj = 1 
is reproduced in the case of no dissipation, 6 = 0. With dissipation, the instability persists 
at high wavenumber. The origin of this instability is uncompensated gas cooling on 6~^ 
crossing times, aided by self-gravitational contraction. In other words, turbulent shocks 
dissipate energy and self-gravity then leads to collapse. 

Figure H] also shows as dashed lines solutions with a dissipation rate proportional to 
6agfik^^'^, which is appropriate for a turbulent gas with a velocity dispersion proportional to 
the square root of the size: Ug oc There is still instability at high k, but the dissipation 

rate for each k is lower when there are slower motions at that k, and so the growth rate is 
lower too. The rest of this paper considers only solutions with ag independent of k, but the 
same shift toward lower growth rate at high k would be present in the full solutions too if 
(Tg; increased with size. 



4. Evaluation of the Dispersion Relation 

4.1. Gas and Stars Considered Separately 

The first term on the right of equation (^^, when set equal to 1, is the dispersion relation 
for pure gas, and the second term is the dispersion relation for pure stars. Combined, they 
are the dispersion relation for stars-|-gas. We show each of these in the next 2 figures. The 
scale heights come from Appendix A. 

Figure [3] shows the pure-gas dispersion relation on the left and the pure-star dispersion 
relation on the right. Both peak at dimensionless wavenumber q ~ 1—2 and both have growth 
rates that increase with decreasing Q. In the case of pure gas, the effect of the dissipation 
term 6qS in equation f l24p is to stretch the range of instabilities to larger wavenumber. When 
6 = 0, the growth rate goes to zero at g ~ 2 for Qg = 0.5 and at g = 3.7 for Qg = 0.3. 

In the case of pure stars (Figure [3] right), the full growth solutions are plotted as solid 
blue curves. These include the scale height term, 1 + qHs with a value of Hs that increases 
from 0.34 at = 0.35 for the top blue curve, to 0.39, 0.49, and 0.58 at = 0.4, 0.5, and 0.6 
for blue curves with decreasing heights. They also include the J-growth term. For comparison, 
the dashed green curve has Qg = 0.5 and no scale height term. The growth is then faster 
because the full stellar column density is used, not the value reduced by 1 + qHg. For the 
value of Hg = 0.49 in this Qs = 0.5 case, 1 + qHs ~ 1.49 at g = 1 near the peak, and this 
makes the growth rate like the one for a full dispersion relation with Qs = 0.5/1.49 = 0.34 at 
the same q. Figure |3] includes this case as the top blue curve (Qs = 0.35), for confirmation; 
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the rest of that blue curve differs from the Qs = 0.5, Hg = case because the scale height 
term also has a q dependence. That is, at increasing q, the term 1 + qHg becomes larger and 
the gravitating effective column density, Ss/(l + fcffs), smaller, making the disk less unstable 
than the Hg = solutions. 

In a second comparison. Figure [3] (right, red dotted line) also shows the dispersion 
relation for Qs = 0.5 and Hg = in the case where the oscillatory solution for the equations 
of motion are used, with J-" from equation (jS]). This solution has no physical meaning because 
at such low Qs the solutions are unstable and not wave-like, but it does show the difference 
that comes from using oscillatory solutions for the reduction factor, rather than growth 
solutions. The oscillatory solutions vary approximately as J-" ~ (1 — sf) /(I — sf + g^), which 
is smaller than the growth solutions which vary as J-growth ~ (1 + ■s^)/(l + + g^). This is 
why the red dotted line in Figure E] is lower than the green dashed line for g > 1. 

Figure H] shows another case with pure gas, this time for Qg > 1, which is normally 
stable. Here the destabilizing effect of the dissipation term, 6qS, can be seen. The solid blue 
curves are a sequence of increasing = 1, 1.2, and 1.4, from left to right, with dissipation 
6 = 1. The dashed blue curve is the same as the Qg = 1 curve but with smaller dissipation, 
6 = 0.5, and so smaller growth rate. The solid red curve has the same Qg = 1 and 5 = 1 as 
the top blue curve, but it has a larger scale height obtained by setting S = OA instead of 
0.2 in the other curves. 

Figure m indicates that although the dissipative disk is unstable for some Qg > 1, the 
growth decreases quickly as Qg increases. Eventually the disk becomes stable at large enough 
Qg if the scale height term is included. Recall in Section |2] we said that infinitesimally thin 
dissipative disks are unstable for any Qg as long as g > Qg/2S. With finite thickness, there 
can be stability at high Qg. Setting s = in equation and using only the first term on 



the right, which is for gas only, we see that if 2gS'/ I Q 



1 + qHg 



> 1, then that equation 



admits unstable solutions with s > 0. This condition rearranges to 

q > 9^^. (28) 

2S Q gH g 

With i^g = we get g > Qg/2S again, but with Hg > 0, we also need Qg < 2S/Hg. 
Unwrapping all the dimensionless variables, this condition is Qg < 2°-^ if if = crg/(7rGIlg^o) 
for a one-component dissipative disk. 
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4.2. Gas and Stars Together 

Solutions to the combined gas+star dispersion relation f l2^ are shown in Figure |5l The 
standard one-component solutions without thickness corrections {Hi = 0) are shown for 
comparison in the left-hand panel as the multicolored dotted line at small q and s. This 
line is made from the superposition in different colors of a pure-star solution with Qs = 0.9 
and Qg = oo, a pure-gas solution with Qs = oo, Qg = 0.9, and 6 = 0, and a gas+star 
solution with gas and stars equivalent, using Qs = 1.8, Qg = 1.8, 6 = 0, and S = 1. They 
are expected to have the same dispersion relation and they do. Also shown in the left-hand 
panel as red curves are two pure-gas solutions {Qs = oo, Qg = 0.9) with different dissipation 
rates, 6 = 0.1 (lower curve) and 6 = 0.2. The extension to high q because of dissipation 
was also shown in Figures [2llll The solid blue curves are two gas+star one-component-like 
cases with Qs = 0.9, Qg = 0.9, 6 = and S = 1, where in the top curve there are no 
thickness corrections {Hi = 0) and in the bottom curve there are. Thickness corrections 
lower the growth rate by about a factor of 60% in this case, where Hs = Hg = 0.44, because 
(1 + qHi)^^=O.Q at q = 1.5 near the peaks of these curves. Finally, another comparison with 
two different dissipation rates is shown by the dashed blue lines but now with gas+stars 
together. Here Qs = 1, Qg = 0.5, and S = 0.5 in both cases, while 6 = for the bottom 
curve and 6 = 0.5 for the top curve. The extension to high q because of dissipation is present. 

A range of values for Qs and Qg are shown for the two-component dispersion relation 
with thickness corrections in the middle panel. The dissipation parameter is fixed at 5 = 0.5 
and the velocity dispersion ratio is 5* = 0.5. These four parameters all enter into the 
dispersion equation fl24|) directly. Also fixed is the ratio of the rotation speed to the stellar 
dispersion v^-ot = 'Omax = 8, the rotation curve slope a = 0.4, and the dark matter to disk 
ratio, D = 1, which enter into the scale height equations flAlip and flA12p . The sequence of 



blue curves with decreasing height has Qs = 0.5, 1, 1.5, and 2 for Qg = 1. The sequence of red 
curves with decreasing height has Qg = 0.5, 1, and 1.5, with Qs = 1. The overlapping case 
has a dashed blue and red curve. These two sequences are slightly different. As Qs increases, 
the peak in the dispersion relation at g ~ 1 — 2 gets less prominent and the dominant part 
of the curve shifts to the right. All that is left to drive the combined instability at low Qs 
is the dissipation in the gas, which extends the instability to high q. On the other hand, as 
Qg increases, the peak at g ~ 1 — 2 remains and even gets relatively stronger as the stellar 
instability begins to dominate. In all cases, there is a tail of unstable growth at high q from 
turbulent dissipation. 

The effect of the dissipation term 6 is shown as a sequence of blue curves in the right- 
hand panel of Figure [5l and the effect of the velocity dispersion ratio S = ag/as is shown as 
a sequence of red curves. For all of these, Qg = 1 and Qg = 1. In addition, -Oj-ot = "i^max = 8, 
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a = 0.4, and D = 1 (which enter into the scale heights and do not matter much). For 
the blue curves, 5 = 0, 0.1, 0.3, 0.5, and 1 with S = 0.5, going from bottom to top, and 
for the red curves S = 0.1, 0.2, 0.3, 0.5, and 1 with 6 = 0.5 going from right to left. The 
overlapping case has a dashed blue and red curve. The effect of increasing dissipation is to 
increase the growth rate at high q, leaving all else about the same. The effect of increasing 
the gas velocity dispersion relative to the stars is to decrease the dimensionless wavenumber, 
i.e., make the unstable wavelength longer. This change is because of the larger gaseous 
scale height as S increases. Along this sequence of S, the dimensionless scale heights are 
{Hg,H,) = (0.052,0.81), (0.10, 0.71), (0.15, 0.65), (0.24, 0.56), and (0.49, 0.49), respectively. 
For the entire sequence of 6 in the blue curves, {Hg, Hg) is the same, (0.24, 0.56). 

Low S also makes the gas+star system stable at low q in the right-hand panel of Figure 
O This is because the gas self-gravity term is small for small S at fixed Qg. Note that the 
ratio qS/ Qg appears in the first term of equation (124|) but S does not appear in the second 
term. As S decreases, the first term becomes smaller compare to the second. For Qg = 1 in 
this Figure [5l the second term alone, from stars, is not enough to drive an instability when 
thickness effects are included. When q is small and S is large, the first term, from gas, adds 
enough forcing to make the system unstable. But when q is small and S is small, the gas does 
not contribute enough. Another way to see this is that the ratio qS/Qg equals ■nG'Lgfik / k'^ 
in physical terms. If S decreases at constant Qs and Qg, then Sg_o/Es,o decreases, so the gas 
becomes less self-gravitating. (Usually one considers that decreasing S corresponds to colder 
gas and that seems like it would be more unstable, but that is only for a constant Sg^o aiid 
not for a constant Qg.) 



4.3. Ratio of Perturbed Surface Densities 

The relative contributions from perturbed stars and gas at the onset of the instability 
may be determined from the ratio of perturbed surface densities using equation fl20l) . which 
involves the quantity Y^. This ratio depends on the wavenumber and growth rate and so 
depends on the solution to the dispersion relation. In dimensionless terms, equals the 
second term on the right of equation (12^ . 

Figure [6] (top) shows Ss/Sg versus dimensionless wavenumber q for all of the gas+star 
dispersion relations shown in Figure El using the same order of panels, curve colors and 
curve types (the previous curves for single component solutions are not plotted in Fig. [6]). 
The bottom panels show the surface density ratio normalized to its initial ratio. Generally, 
the stellar contribution to the instability is proportional to the initial star-to-gas ratio. It is 
equal to this ratio for g < 1, and less than it for g > 1. 
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For the fastest growing wavenumber, which is usually larger than g ~ 1 (Fig. [5]), the ratio 
of stars to gas in the perturbation can be between 0.5 and 1 times the initial ratio. Figure 
[7] shows histograms of the normalized ratio of perturbed stars to gas at the wavenumber of 
peak growth. Each histogram is taken from a series of 800 runs with Qs = 0.2, 0.4, 0.6, up to 
4, and Qg = 0.1, 0.2, 0.2, up to 4, and for the cases in this range where there is an instability 
and where the dimensionless wavenumber of peak growth is less than 40. The dissipation 
parameter S increases from top to bottom in the figure. With more dissipation, the peaks in 
the histograms of the left-hand panels, which are for S* = 1, shift to the left, which means 
there is an increasing extension of SgSg^o/SgSs^o to very low values. These changes mean that 
the contribution from gas is becoming more important as turbulent dissipation increases. 

The right-hand side of Figure [7] is for lower relative gas velocity dispersion, S = 0.2. 
Aside from an increase in nearly pure-gas instabilities at SsSg o/SgSs^o ~ with increasing 
6, there is not much change in the position of the peak. The lower cutoff to the main peak 
at SsSg^o/SgSs,o ~ 0.4 means that gas contributes less to the instability when S is low than 
when S is high, even when there is dissipation. The reason for this is that when S is low, 
the first term on the right in equation fl24p is low at fixed Qg. 

To derive the absolute ratio of perturbed stars to gas, we recall that the unperturbed 
ratio is Ss,o/Sg^o = Qg/SQs- Then the absolute ratio is Qg/SQs times the relative ratio in 
Figure [71 For S = 0.2, Figure [7] has a peak at a relative ratio of ~ 0.6, and so the absolute 
ratio of perturbation strengths is 

|i = §ix^~3|. (29) 

For S = 1 and 6 = 0, Figure [7] has a peak near 0.9, so Ss/Sg =~ 0.9Qg/Qs. For S = 1 and 
6 = 1, the peak is near 0.6, so Sg/Sg ~ 0.6Qg/Qs. Generally Qg < Qs, so the perturbed 
mass in stars is less than or comparable to the perturbed mass in gas for a wide range of 
cases. When the Q values are large and the main instability involves strong gas dissipation, 
the normalized density ratios in Figure [7] are much lower, ~ — 0.2, and then the relative 
contribution from stars is less. 

A contribution from stars that is comparable to the contribution from gas in an insta- 
bility, Ss/Sg ~ 1, implies that spiral arms which form from the instability contain old disk 
stars. This is obvious and generally acknowledged for spiral arms. However, the same high 
fraction of disk stars could also be present in instabilities at the dominant wavenumber that 
make large globular clouds, as in dwarf galaxies with l ow shear or in highly unstable clumpy 



galaxies in the early universe flElmegreen et al.ll2009[ ). The oldest disk stars would be less 



abundant in these clouds than the youngest disk stars because the stellar velocity dispersion 
increases with age. Still, stars of intermediate and young age that did not form in the clouds 
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could have an excess abundance compared to the surrounding disk. According to equation 
f l29|) . this excess in disk stars approximately equals the excess in gas multiplied by the ratio 
of Qg to Qs. 



4.4. Instability Thresholds 



The stability threshold is given by the dispersion relation at a growth rate of s = 0. 
In the absence of dissipation {6 = 0) and thickness effects {H = 0), equation flMj) is a sum 
of two terms, 2qS/ (Qg [1 + q'^S^]) and 2q/ {Qs [1 + q"^]), considering that equation f[T5]) is a 
good a pproximation . Each term has a maximum over variations in q equal to 1/Q. This 
is how I Wang fc Silkl (119941 ) derived an approximate instability threshold for the two-fluid 
instability: 



1 



> 1. 



(30) 



<5tot,ws Qg Qs 

This threshold has been widely used and is more conve nient than the exa ct Q threshold for 
these conditions that requires solving a cubic equation (lElmegreenl Il995l ) . Extrapolation of 
this method with thickness corrections might suggest that instability of a more realistic thick 
disk requires 

11 1 

(31) 



1 



+ 



1 



Qtot.wsH Qg{l + qHg) Qs{l + qH, 



> 1. 



Recently, iRomeo fc WiegertI (120111 ) proposed a modified form of the Wang & Silk equa- 
tion that is more accurate in the standard case (^ = 0), 



1 W 

T^T + 7=r > 1 if Qs> Qg 



Qtot.RW Qg Qs 

w 1 
= 7^ + 7^ > 1 if Qs < Qg. 



The weighting function is 



W{S) 



2S 



(32) 
(33) 

(34) 



1 + ^2' 

where S is the ratio of velocity dispersions, cxg/as. The weighting function allows the gas 
term in equation (12^ to contain S. With thickness effects. 



W 



> 1 if Qs> Qg 



Qtot,RWH Qg{l + qHg) Qs{l + qHs 

W 1 
= \ — > 1 if Qs < Qs- 

Qg{l + qH,) Qs{l + qHs) 



(35) 
(36) 
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We recall here for convenience that the dimensionless quantity qH is identical to the physical 
quantity kH. We test these approximations by finding their values at the limit where s = 
in equation We also test how well they work when gas dissipation is included. 

Figure E] shows numerous plots of peak growth rate versus l/Qtot- Each panel contains 
20 blue curves and 20 red curves. The 20 curves correspond to values of Qs = 0.2, 0.4, ...4, 
and, within each curve, Qg varies from 0.1 to 4 in steps of 0.1. Thus each curve is a sequence 
of decreasing peak growth rates for increasing Qg. The peak growth rate is evaluated only 
up to a wavenumber of g = 40. 

The left-hand panels plot the Wang & Silk parameter on the abscissa, with the blue 
curves for equation fl30|) and the red curves for equation fl3T]) . The middle column plots the 
Romeo & Wiegert parameter on the abscissa, again using blue to indicate no thickness terms 
in that parameter. These two columns use a ratio of dispersions S = 0.2, which is typical for 
spiral galaxies. The right-hand panels consider 5* = 1, for which the Wang & Silk parameter 
is the same as the Romeo & Wiegert parameter. 

Thickness effects are considered in the dispersion relation solutions for all of the curves 
in the top four rows (even when these effects are ignored in the expression for 1/Qtot)- What 
differs in these top four rows is the value of the dissipation parameter, 6. The bottom row 
solves the dispersion relation without dissipation and without thickness terms, taking 
Hs = Hg = 0; the red and blue curves are the same then. 

In all panels, the threshold of stability occurs where the curves intersect the x-axis, at 
Speak = 0. Evidently, the stability thresholds given by equations (!30|) - (!33|) work well in the 
situations for which they were designed. For an infinitesimally thin disk with no dissipation 
and a velocity dispersion ratio of S* = 1 (bottom right panel), the Wang & Silk formula 
works perfectly: all the curves intercept the l/Qtot,ws = 1 point at Speak = 0. In this limit, 
the two fiuids act as a single fiuid because they have the same velocity dispersion. When 
5* = 0.2 at zero thickness and no dissipation (lower left panel), the stability threshold is 
l/Qtot,ws = 1-28 on average for all intercepts. The Romeo & Wiegert relation works well for 
both S values in the 5 = H = limits (i.e., in the bottom middle and bottom right panels). 

When thickness effects are considered in equation fl24|) . but still without dissipation (top 
panels), both the Wang & Silk parameter (red curves: top left and top right panels) and the 
Romeo & Wiegert parameter (red curves: top middle and top right panels) work well in the 
sense that they intercept the x-axis at 1/Qtot ~ 1- Thus equations ( ET\i and (15^ are both 
reasonable approximations to the stability limit on average. 

Dissipation lowers 1/Qtot in all cases. Both the red and blue curves in Figure [8] shift 
to the left in a sequence of increasing S going from the top row to the fourth row. With 
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zero thickness, disks are always unstable for large enough wavenumber when dissipation is 
included (see discussion after equation |2l]), but for finite thickness, as in the top four panels 
of Figure El a disk eventually become stable when 1/Qtot is small enough (see discussion after 
equation The blue curves dive into the abscissa at low l/Qtot,ws and Qtot.RW because of 
this disk stability. The red curves approach the abscissa asymptotically as l/Qtot,wsH and 
l/Qtot,RWH go to zero because q at peak growth gets very large when there is dissipation. 
Then (1 + qH) pulls down the l/Qtot parameters at large Q. 

The irregular excursions in the red curves of Figure [Dare because of a rapidly decreasing 
wavenumber at peak growth for increasing Q at very low Q. Recall that for a one-component 
instabihty, /cpeak = irGTi/a"^, or gpeak = {QS)~^. This large q at small Q enters into the 
(1 + qH) factor, and, this, along with slight variations in H, makes the excursions. 

Figure [3 shows the average x-axis values where the growth rates in Figure [S] are between 
0.03 and 0.3. This is a reasonable range to indicate either absolute stability or slow growth in 
the cases of asymptotic approach to the abscissa. Blue curves are again for Qtot parameters 
that ignore thickness and red are for those that include thickness (where all of the dispersion 
relation solutions include thickness). The open square at 5 = is the Wang & Silk equation 
fl30|l for 5* = 0.2 and i^^ = in the dispersion relation; the x-symbol is the Romeo & Wiegert 
result for 5* = 0.2 and H = 0, and the open triangle is for both Wang & Silk and Romeo & 
Wiegert at 5'=1 with H = 0. Figure |9] confirms the result 

l/Qtot parameters decreases with increasing dissipation: for 6 in the range from 0.4 — 1, 
Qtot,wsH and Qtot,RWH should be less than ~ 1.5 — 2 for instability when S = 1, and they 
should be less than ~ 3 — 4 for instability when S = 0.2. If thickness effects are ignored in the 
stability condition (blue curves), then for this same range of 6, Qtot.ws = Qtot,RW < 0.78—0.84 
for instability when S = 1, and Qtot,ws < 0.70 — 0.74 and Qtot.RW < 0.89 — 0.95 for instability 
when 5* = 0.2. 

These results indicate that unknown dissipation partly compensates for unknown thick- 
ness in the stability condition. That is, Qtot conditions that are evaluated incorrectly without 
including the thickness terms (blue curves in Figure [H]) have a threshold for instability that is 
still close to 1 when dissipation is important. This compensation has been evident for a long 
time in simpler derivations where one uses the adiabatic expression for pressure, P oc p'^, 
instead of t he perturbed energy equation ( ITTj) . Denser regions are cooler so 7 < 1. The com- 



pilation in iMyerd (119781 ) suggests that 7 ~ 0.25 at interstellar densities between 0.1 cm~^ 
and 100 cm~^. This index appears in Qg as 7^/^ multiplying the velocity dispersion, cXg. Thus 
a version of Qg with adiabatic cooling and thickness corrections is 7^/^(TgK(l + kH)/7iGT,g. 
Since kH ~ 1 usually, 7^/^(l + A;if ) ~ 1, indicating that cooling effects offset thickness effects 
even in a simple derivation. 
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5. Azimuthal Instabilities 

The main result of the previous section is that dissipation promotes gravitational insta- 
bilities at high wavenumbers even when a gas+star disk is stable by conventional criteria. 
Thickness effects ultimately stabilize the disk for high Qg and Qs- In the mos t extreme re- 
ductio n of this result to a single stability parameter Qtot, we find that both the lWang fc Silk 
( I1994J ) and the lRomeo fc Wiegertl ( 1201 ll ) definitions are reasonably good provided the thresh- 
old is raised by a factor of 2 to 3, depending on the relative dissipation rate. For example, 
if significant energy dissipation in the gas occurs in 2 crossing times over a perturbed wave- 
length (5 = 0.5), then Qtot including thickness effects has to exceed about 2 for stability, 
instead of the usual value of 1. This is for radial instabilities. 



A higher threshold is also likely for azimuthal instabilities of the swing-amplifier type 
(iToomrdllOSlf). There is a second instability parame ter in that cas e ., i.e. , X in the notation 
of I Julian &: Toomrd (Il966[ ) and J in the notation of iLau fc BertinI ( 119781 ). 



The X parameter is defined as 



X 



At 



A 



(37) 



crit 



where At is the wavelength of the perturbation in the azimuthal direction and Acrit = 27r//ccrit 
for fccrit = K^/(27rGS). The wavelength in the azimuthal direction is 27t R/m for radius R 
and number m of perturbations in the azimuthal direction. Thus 



X 



Rk' 



RQ Q 



27TGT,m 2m 2q^ 



(38) 



where quantities have been normalized in the same way as before: R = Rn/a and qt = kta/n, 
where kt = 2-n /\. 



The J parameter is about the same: 



27iGj:\ f2n\ 



^2 ; 



dlnQ 
dlnR 



1/2 



1 /2n\ 
x\~) 



dlnQ 
dlnR 



1/2 



(39) 



The terms in parentheses on the right are of order unity for angular rotation rate Q[R). For 
a fiat rotation curve, their product is 1.4. 



Toomrd (119641 ) noted that the azimuthal instability in a rotating galaxy has the same size 
limit as the radial instability, i.e.. At < Acrit, which is the requirement that gravity overcomes 
the centrifugal force that arises when condensing material spins up by the Coriolis force. This 
means that instability is favored by small X or large J. 
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Lau fc BertinI (119781 ) and iBertin et al.l (Il989[ ) introduced an effective two-dimensional 



stability parameter 

1 1 



'eff 

The second term is 

72^4 /on\ 2 



4 U^IJ-''"®' 

where the factor 2 is for a flat rotation curve. 

The azimuthal instability is stronger if the azimuthal wavenumber gt is larger. There 
is no absolute stability threshold if arbitrarily large gt (or "^) are allowed, because Qeff can 
always be made less than the threshold value of 1 by large enough q^. However, at high 
there is also pressure stabilization for an adiabatic gas undergoing azimuthal instabilities. 
This is the usual Jeans condition. Without an adiabatic gas, this Jeans condition goes away 
if turbulent dissipation is fast enough, as shown in previous sections. 

Transient growth of shearing spiral waves is still possible. If the dissipation rate is 
faster than the shear rate, then transient growth can lead to runaway growth. What usually 
stabilizes the swing amplifier at late times is high pressure at the high wavenumbers that 
are reached by the swept-back, tightly-wound wavelets. With turbulent dissipation, this gas 
pressure becomes small. Stabilization at late times then comes from thickness effects rather 
than pressure. That is, thickness dilutes the gravitational force at high wavenumbers by the 
factor (1 + kE)-^. 



Discussion 



6.1. Pervasive Instabilities 



Turbulent dissipation, viscosity, and magnetic fields destabihze galaxy disks. All three 
either decrease or remove the threshold column density for instability. Dissipation does this 
by removing the stabilization from turbulent pressure, as shown here, wh ile magneti c fields 
and viscosity r e move the stabilization from Coriolis f orces, as s hown by lElmegreenI (119871 ) 
and iKim et al.l (120021 ) for the magnetic case, and by iGammid (119961 ) for the viscous case. 
All three processes rely on gas. The effects of dissipation discussed here arise if the gas is 
not isothermal or simply adiabatic, but satisfies an energy equation with excess dissipation 
over heating in a perturbed region. The particular form for the dissipation rate assumed 
here, i.e. proportional to the initial crossing rate in a per turbation, ba^k, \s, appropriate for 



supersonic turbulence in the general interstellar medium (IMac Low et al.lll998l : IStone et al. 
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19981 ). This result differs qualitatively from the usual result that cooling makes a gas disk 
more unstable by lowering Qg. It also differs from the result that adiabatic gas corresponds 
to greater instability when the adiabatic index 7 is lower. 

Our primary concern is the stabil i ty of g alactic gas. This s ituation differs fro m the 
rotating gas disks studied by iGammid (I2OOII ). iRice et al.l (120031 ) . iMejia et al.l (120051 ) . and 
others in which gravitational fragmentation occurs if the dissipation time is less than about 
half the rotation time. These authors assume a constant cooling time Tc for an energy loss 
rate U/ with internal energy density U. This would take the place of our term F — A in 
equation (|9]) and therefore differ from our assumed rate, 5agkP/{^ ~ 1); in equation (|T0|) . 
In our equations, the perturbed cooling time is not constant but scales with the size of the 
region, oc 1/k. The other studies consider instabilities in optically thick, subsonic accretion 
disks around protostars or active galactic nuclei. In this environment, dissipation is difficult 
because the radiation c annot easily escape. Dissipation may be aided by convection in the 
l ower layers (IBossI 120041 ). or it may be prevented by slow radiative losses at the disk surface 
( ICai et al.ll2010l ). Galactic gas, in contrast, is optically thin and supersonic. Dissipation is 
by random shocks and energy loss is by unblocked radiation from the shock fronts. In our 
case, small-scale instabilities are possible in the rapidly cooled shocked regions, even if they 
are smaller than the Jeans length at the initial velocity dispersion. 

The present result implies that galaxy disks may be more unstable than previously 
determined from models or computer simulations with isothermal or adiabatic gas. The 
effective Q threshold for gas+star radial instabilities increases to ~ 2 if turbulen t dissipation 
takes a few crossing times. This increased Qtot value is comparable to what iLeroy et al. 



( I2OO8I ) found for most regions of their galaxies, which implies that galactic disks are only 
marginally stable. 

Spiral arms are ex amples of gravitational in stabilities that depend on Qtot- Multiple 
arm galaxies have what iToomre fc Kalnajd (Il99l[ ) called spiral chaos, with gas+star spiral 
features of all sizes forming everywhere in the disk and continuously replacing old ones that 
get sheared away. The gaseous parts of disks can have inde pendent instabilit i es th at make 
spiral arms in a fiocculent pattern of star formation (e.g., lElmegreen et al.ll201ll ). Mul- 
tiple arm and fiocculent structure are instabiliti es covering a wide ran ge of scales. They 



should pump turbulent energy into the ga s (e.g., iBournaud et al.ll2010l ) and heat the stars 



(ICarlberg fc Sel 



wood 



1985 



Sellwood 



20111). Grand design spirals may be triggered by inter- 



actions or bars (IKormendy fc Normanlll979l ) but require a marginally stable disk to amplify 



these perturbations to large arm/interarm contrasts. 

Power spectra of optical light in galaxies, including grand-design galaxies like M81, 
are continuous power-laws suggesting pervasive instabilities and a cascade of power over a 
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wide range of scales ( lElmegreen et al.ll2003l 120061 ) . There is no dominance of structure at 
the fastest-growing wavelength of an instability. Such lengths are a feature of the linear 
regime and may appear only in young regions that start fairly smooth. Examples might be 
expanding shells or spiral arm shocks, which sometimes show a characteristic clump scale 
(iDeharveng et al.l l2010l : lElmegreen et al.l |2006| ) . 



6.2. Feedback Regulation of Disk Stability 



The observation by iLeroy et al.l ( 120081 ) that Qtot ~ 2 for most galaxy disks in their 
survey implies some regulatory mechanism. Regulation is usually considered to involve 
star formation so that low Qt.ni, triggers more star formation, which increases cTg and Qtot 
( iGoldreich fc Lynden-Belll Il965l ). High Qtot triggers less star formation, causing cTg to de- 
crease by passive cooling, and lowering Qtot to a more unstable value. Qtot also affects the 
stellar velocity dispersion because growing instabilities scatter disk stars. The present results 
detach cXg from Qtot somewhat, because dissipation also lowers during the initial growth 
of the instability. Dissipation and instabilities are not two separate steps in this model. 

The variety of spiral arm strengths in galaxies suggests that Qtot is not strongly reg- 
ulated. The stellar spirals in multiple arm and fiocculent galaxies vary in strength, from 
50% or stronger surface density modulations in multiple arm galaxies to no perceptible 

Of the 197 fiocculent galaxies catalogued by 



Elmeereen & Elmeereen 


( 


1987 


). 85% 


2MASS survey ( 


Elmegreen et al. 


2011) 



more stable stellar disks than multiple arm galaxies, which means that not all Qtot values are 
the same. Similarly, in a sample of 147 Ohio State University Bright Galaxy Survey galaxies 
less inclined than 65°, the distribution functi on of spiral arm torques is broad and featureless, 
spanning a factor of ~ 4 (iButa et al.ll2005l ). Spiral arm morphology and strength are the 
most direct indicators of Qtot, because these structures are what gas+star instabilities make. 
Star formation does not have such a clear connection with Q 



tot- 



The composite stability conditions shown in Figure [8] illustrate the conventional notion 
that instabilities maintain Qtot at a marginally unstable value by stirring the stars and gas. 
If Qtot is small then the threshold, Speak, is large and stirring is active, which lowers Qtot- 
Feedback is complex for a gas+star fluid, however, because the stars do not cool. If star 
formation is active, then cool stars can be ad ded continuously to the stellar population, 
which is like cooling (ICarlberg fc Sellwoodlll985l ). but the same instabilities also heat the old 
stars in proportion. Stronger instabilities heat more old stars and, through cloud formation, 
produce more cool stars. If heating beats cooling in this proportion, then there would be no 
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feedback to make the stars cooler if Qtot is too high. If coohng from new stars beats heating, 
then there would be no feedback to make the stars heat up if Qtot is too low. 

Feedback stabilization of Qtot is also complex for a gas+star fluid because the cool young 
stars have the same velocity dispersion and mass as the gas did before the stars formed. Thus 
the new, cool stars do not add any extra instability to the combined disk. They remove the 
ability to dissipate energy however, which is a stabilizing effect. 

Persistent heating of old stars by gravitational instabilities requires more and more 
cool gas to maintain a marginal Qtot if there is regulation. This eventually leads to a gas 
instability that is decoupled from the star instability, i.e., when S = cXg/as << 1. Then there 
are two peaks in the dispersion rela tion, one for stars at low q and another for gas at high q 



(jjog fc Solomon! 1 1984J : lRomedll992l ). Some flocculent galaxies could be in this regime. 



If the stellar Qs is low, then stellar spirals form readily and the disk cannot be stabilized 
by extra turbulence following star formation. Feedback comes only from stellar heating in 
the waves. On the other hand, if Qs is high, then the combined instability shows up primarily 
in the gas, producing flocculent spirals and filaments of star formation. 

Galactic gas accretion can modify Qtot by adding cold gas, and after star formation, by 
adding slow-moving stars. Following a sudden influx of gas, Qtot should drop and a period 
of intense spiral arm or clump formation and star formation should begin. When the gas 
eventually gets converted to stars, Qtot should increase and the spiral arms weaken. This 
would be a regulatory process that does not connect Qtot with star formation directly but 
still keeps Qtot at a more-or-less constant value, particularly if the star formation rate is 
about equal to the galactic accretion rate on average. 



6.3. Star Formation 



6.3.1. Is there a Connection with Galactic-Scale Gravitational Instabilities? 



The dir ect role of g r avita ti onal instabil i ties i n star formation is difficult to observe. 
According to lLeroy et al.l ( 120081 ). iBigiel et al.l ( 120081 ) and others, star formation only requires 
molecular gas, and molecular gas only requires sufficient pressure, most of which comes from 
the weight of the gas layer in the disk. There is nothing about cloud formation or possible 
instabilities that induce cloud formation i n these conditio n s. On ly pr essure determines the 
abundance of dense, star-forming clouds. iKrumholz et al.l ( l2009l ) and lOstriker et al.l (120101 ) 
reproduce the observed radial profiles of galaxies using only the constraint that molecular 
clouds are shielded by atomic gas in some type of equihbrium. It is difficult to see how 
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galactic-scale gravitational instabilities can influence these results as the instabilities operate 
on the total gas (and stars) without regard to molecular phase. If gravitational instabilities 
directly triggered star formation, then the primary scaling law should have a star formation 
rate that depends on a high power, such as 1.5 or 2, of the total gas column d ensity (e.g.. 



Madore Il977 



2007 



Li et a. 



LiuetaL 



20051 ). and this is apparently not observed fhowe ver. see lYang et al. 



(]2007f ). 



20111 ). For a review of star formation processes, see iMcKee fc Ostriker 



Perhaps this disconnection between star formation and gravitational instabilities makes 
sense in retrospect. If the gas is already highly molecular, then the conversion of a little extra 
atomic material into molecular material by compression in instabilities will not influence the 
total molecular mass and star formation rate much, nor will any rearrangement of the existing 
molecular clouds by the unstable motions. Even cloud collisions during an instability would 
not a ffect the star formation rate if the molecular depletion time is as constant as it appears 
to be iLeroy et al.l ( l2008l ). On the other hand, if the gas is highly atomic, then instabilities 
could promote molecule formation if the density contrast becomes large. However, on the 
dominant scale of the gas+star instability, which is (J^/ k (i.e., q ~ 1), and for the dominant 
growth rate, ~ k (i.e., s ~ 1), these instabilities would only increase the density by several 
exponential growth factors during the available shear time. This is not enough compression 
to make molecular clouds. Star formation has to occur in the non-linear regime of unstable 
growth, or independently of this growth, when the ga s separates from the st ars and collapses 



into dense molecular clouds by its own gravity (e.g., iKim &: Ostrikerll2002l ). The stars that 
helped start the instability will have moved away in their epicycles after this time, with little 
long-term excess following the clumped gas. Star formation also occurs after the instability 
causes contraction in the third dimension, considering that giant molecular clouds (GMCs) 
are smaller than the scale height. In addition, star formation can occur in pressurized shells 
and cometary structures that are unrelated to galactic-scale instabilities. 

The observation of a linear molecular-cloud star formation law, which suggests a sec- 
ondary role for galactic-scale gravitational instabilities, combined with the observation of a 
weakly- varying Qtot, which suggests feedback-regulation of Qtot by gravitational instabilities 
that may trigger star formation, poses an interesting contradiction. Evidently the role of the 
these instabilities in star formation is indirect. 

One possibility is that gravitational instabilities replenish the self-gravitating molecular 
gas that gets converted into diffuse gas during the star formation process. This way the 
molecular gas abundance remains constant on average, and the star formation rate does not 
scale with the cloud formation rate, but only with the current mol ecular abundarice. A n 
example is provided by the simulations in iBournaud et al.l (120101 ) and iHopkins et al.l ( 1201 ll ) , 
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which show gravitational instabilities that form self-gravitating gas from diffuse gas, and 
star formation that disrupts these gravitating clouds, converting a fraction of their gas back 
into diffuse form. The conversion rate of diffuse gas into self-gravitating gas by galactic-scale 
instabilities can be much larger than the conversion rate of self-gravitating gas into stars by 
local collapse. A simple model for this is discussed in the next sub-section. 



6.3.2. A Simple Model of Star Formation with a Role for Gravitational Instabilities 

A simple model of star formation illustrates how gravitational instabilities can play an 
important role in cloud formation, which is cloud re-assembly in this picture, and yet not 
trigger star formation directly or have a direct connection to the star formation rate. In this 
model, there arc three phases of interstellar matter: self-gravitating clouds in which stars 
form and which are generally molecular, diffuse molecular clouds in which stars do not form 
because their density is too low, and diffuse atomic or ionized clouds, which are also too low 
in density to form stars. If, in a certain region their masses are 7713, m2, and mi, respectively, 
then we can write the equations governing their changes: 

dm^ mi -|- m2 d^m^ 



dt Tqi Tburst Tc 



(42) 



consume 



dm2 _ _ m2 _^ d2m^ ^^^^ 
dt Tgi 

dmi _ _rni_ ^ dims ^^^^ 

dt Tqi Iburst 

The timescales and terms in these equations are as follows. Tgi = 1/^ is the timescale for 
gravitational instabilities studied here. For most cases, this is comparable to the epicyclic 
time, which means that s ~ 1 in the dimensionless notation used above. Then the first 
equation says that the self-gravitating mass increases on the timescale of the gravitational 
instability with the mass coming from diffuse gas. That is, gravitational instabilities com- 
press all of the gas into discrete self-gravitating clouds, but the practical effect of this for 
star formation is to move the diffuse gas that lies between existing GMCs onto those GMCs, 
and to make new GMCs. 

^burst is the duration of a single event of star formation in any one cloud, such as the 
lifetime of an active OB association. It should be something like 10-20 Myr for a 10^ Mq 
cloud. The factor ^3 is the fraction of a self-gravitating cloud that is converted into diffuse 
gas during each burst; it might be 30% or more, considering disruption of a cluster-forming 
core and all of the associated ionization. It is not likely to be 100% because dense clouds 
are mostly pushed to the side by pressures from ionization and stellar winds, into shell-like 



- 25 - 



and cometary structures, without complete disassociation into atoms. Of this disruption 
fraction, part of it presumably stays in molecular form but at too low a density to form stars 
immediately. It could be molecular without r auch CO, for exampl e. Many local diffuse 



clouds are in this state (ISpitzer fc Jenkins! Il 9751 ). iRahman et al.l ( 120101 ) discuss observations 



of the star formation law with explicit consideration of diffuse molecular gas. It could also 
be in the form of GMC envelopes that are exposed to background starlight and have high 
ionization fractions and long magnetic diffusion times (Elmegreen 2007). With this term 
in equation fH21) . the self-gravitating mass decreases because of cloud disruption on the 
timescale of the individual burst. 

^'consume IS the cousumptiou time for self-gravitating gas in general. The consumption 
time is the total self-gravitating molecular cloud mass in a galaxy divided by t he total star 



formation rate. It is u sually measured to be ~ 2 Gyr for main galaxy disks (ILeroy et al. 



20081 : iBigiel et al.l 120081 ). The star formation rate is this last term, m^/T^ 



consume • 



There are also parameters in these equations that indicate the fractions of the various 
phases that convert into other phases. The parameter d2 is the fraction of the gravitating 
cloud mass that turns into diffuse molecular gas during a burst, and di is the fraction that 
turns into diffuse atomic or ionized gas; then di + 62 = d^. 

Equations fj42l) - fH4l) can be solved numerically for the time dependence of the three 
phases, but the result of interest here comes more easily by considering the dominant terms. 
In equation fH21) . the last term represents the conversion of gas into stars. Because of the low 
efficiency of star formation, which means the large value of Tconsume compared to Tburst and 
Tgi, this third term is relatively small. In that case, there is a quasi-equilibrium between 
the three phases in which the time derivatives are small compared to each separate term 
involving Tqi and Tburst- In this equihbrium, 

mi ~ — di ; m2 ~ — d2. (45) 

^ burst ^ burst 

What we would like is a measure of how active the gravitational instabilities are in rebuilding 
the clouds that form stars. The building rate is (mi +7712) /Tqi, because the instability collects 
all of the gas and places it in dense clouds, which means it converts the diffuse gas into self- 
gravitating gas. This building rate should be compared with the star formation rate, which 
is ma/Tconsume- The ratio is 

mi + m2 \ / ^consume \ d-^Tci 



' consume 



Tqi J \ ^3 J Tburst 

For d^ ~ 0.3, Tconsume ~ 1 Gyr, and Tburst ~ 20 Myr, the ratio is 15. 



(46) 



The point of this simple model is to show how star formation can happen in self- 
gravitating clouds at a fixed rate, independent of gravitational instabilities, i.e., without 
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direct triggering, and yet these instabilities are important in maintaining the cloud popula- 
tion so that this star formation processes can continue. The mass flux into self-gravitating 
clouds is 15 times the mass flux into stars in this example. In practice, other processes 
of dense cloud formation, such as shell formation around superbubbles, will combine with 
gravitational instabilities to return diffuse debris back into self-gravitating clouds. 

There are many computer simulations of star formation in galaxy disks that contain the 
basic elements discussed here, namely gravitational instabilities that drive turbulence and 
dense cloud formation, star formation in the d ense clouds with a certain fixed or density- 
dependent rate, and dense clo u d disruption (e.g 
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2009 





Dobbs fc Pringle 



2009 



Baba et al 



2009 



Gnedin &: Kravtsov 



Wada fc NormaJ2007l: iRobertson k. Kravtsov! 



2010 



Bournaud et al. 



Murante et al 



Tasker 



2010 



Gnedin &: Kravtsov 



2011 



Ostriker fc Shettvll201ll : iTasker fc Tan 



201 ll ). These generally have an equation of state for the gas that includes the 



heightened instability from dissipation that is discussed in the present paper. 



7. Summary 

Gravitational instabilities in a gas+star galaxy disk are enhanced when gas dissipation is 
allowed to occur during the growth of the instability. Then the restoring force from pressure 
and the velocity dispersion are reduced during the unstable growth. This has the effect of 
removing the minimum wavelength (the Jeans length) from the problem, and for a thin disk, 
promoting an unstable condition regardless of Q. For a thick disk, the thickness itself is what 
limits or prevents instabilities by diluting self-gravity in the in-plane direction compared to 
the Coriolis force. 

This paper derived the dispersion relation for small perturbations in a dissipative gas + 
stellar disk. The effects of disk thickness were considered in detail, including a self-consistent 
evaluation of the thicknesses of the two components. Threshold effects for an appropriately 
defined parameter Qtot were determined and found to differ from the usual results determined 
without explicit inclusion of dissipation. That is, Qtot has to be a factor of 2 to 3 higher than 
the usual value of unity before a reasonably thick disk becomes significantly stable. This 
change implies that real galaxy disks are more unstable than previously believed. 

The greater instability has implications for spiral waves and cloud formation, and for 
feedback control of Qtot, which is much more complicated for a two-component disk than 
a one-component disk, considering that an existing stellar component can never cool. The 
instability also has implications for star formation laws, particularly if they suggest a linear 
dependence between the star formation rate and the molecular abundance. Such laws would 
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not seem to have a role for gravitational instabilities in directly triggering star formation, 
but they may still have a role in re-assembling the self-gravitating clouds that are torn apart 
by stellar feedback. The mass exchange in this assembly phase can exceed the mass flow into 
stars by a factor of 10 or more, indicating the importance of the instability as a component 
of the whole process. 
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A. Scale Heights 



The vertical profiles for perturbed gas and stars in our model are the same as the vertical 
profiles for equilibrium gas and stars because all of the motions are assumed to be parallel 



to the dis k plane. The se a. 



e hei ghts then come from the equations of vertical equilibrium. 



We follow iNarayan fc Joel (120021 ) and write the equation for density of each component i 
g and s: 

d'^Pi P-^ ( A ni , ^ , dK^M \ , 1 /rfpi 
^ = ^ -4-G(p, + p.) + ^ + - - 



(Al) 



This preprint was prepared with the A AS lAT^rpC macros v5.2. 
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where the dark matter acceleration per unit length is 
dK^M vL^ z'^R^^ 



max 



1 h 



Rc 






r 







arctan [r/i?c] (A2) 



for a dark matter density distribution 

1 /2 

with central density pdmo, core radius -Rc, 3D radius r = {R'^ + z^) , disk radius in cylin- 
drical coordinates i?, and extrapolated rotation speed at infinity, fmax- 

The rotation curve enters into our primary variable k, which is used for normalization, 
so we should be self-consistent in the definitions of these quantities. We write the rotation 
curve in the midplane as the sum of contributions from the dark matter and the disk, 

^rot(^) = <tM^) + ^rot,Disk(^) (A4) 

with 

R 
and 

M^m{R)= I 47ri?VDMC^^ = 47rpDM0^n ^ - arctan— . (^6) 



L,DM(fi) = (A5) 



We assume that the disk mass column density profiles are exponential, Si^o = Si^ooe 
and have the same disk scale lengths Rd for gas and stars. Then 



D 

5 



-'rot 



Di3k(^) = 4vrG (Sg,o + S3,o) RW {h [v] K, [y] - h [y] [y]) (AT) 



in the notation of Binney & Tremaine (2008, eqn. 2.165), with modified Bessel functions / 
and and y = R/2Rd. 

Note that as i? — )■ oo, v^^^ — )■ AttGpdmoRc, which is v^^^ in equation (lA3p . Now we write 

\ R J \ dlogR J ^ ' 

and define a = dlogVrot/dlogR. This gives a relation between the dimensionless disk radius 
R and the dimensionless rotation speed, 

R' = 2{l + a)v^,,. (A9) 
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The other dimensionless parameters are related also. From equation flA3|) . the ratio of the 
dark matter central column density to the stellar disk central column density is 7rpDMo-Rc/Ss,oo- 
If we let D = Pdmo-Rc/Ss,oo, then 

Rc = (AlO) 

In terms of the characteristic dimensions used for section [21 namely k for inverse time 
and o"s for velocity, the vertical equilibrium equations in dimensionless form become 

d'^Pg _ Pg ( 2 dK^M \ 1 / dpg ^ ^ 



h(-l(h^P.)^'-^]^l('4Y- (A12) 



dz'^ \ Qs dz I Ps \ dz 

Equation (]A2|) for dK^u/dz looks essentially the same in dimensionless form with dimension- 
less distances like f = rn/os substituted for physical distances, and with the dimensionless 
rotation speed v^ot = '^rot/o's- In equations (lAlip and (lA12p . dimensionless density is related 
to physical density as p = 2ttGpQs/ n"^- Note that for a single fluid, the solution to equa- 
tion flAllI) would be p(z) = p(z = 0)sech^ (z/H) and the dimensionless scale height would 



z 



he H = S (Qs/ p[z = 0]) . The dimensionless column density would be So = j°^ooP^ 
2p{z = 0)H. For a single fluid in physical units, H = a^/ (vrGSo) = cx/ {2'KGp[z = 0]Y^^. 

The vertical equilibrium equations may be solved for Pi{z) and the scale heights defined 
as Hi = 0.5Sio/pi(;z = 0). For the dimensionless solutions, the dimensionless input param- 
eters are Qs, S = dg/cis, -Orot, Vma.x, and the dark matter-to-disk column density ratio, D. 
Also entering is the logarithmic gradient of the rotation curve, a, which varies with radius 
and depends on D, and disk-to-halo size ratio, Rq/Rc- The unperturbed gas-to-star column 
density ratio, Sg o/Sg^o niay be treated either as an independent variable, or it can be taken 
equal to SQs/Qg if Qg is independently specified. 

In practice, we do not need a detailed model of a galaxy just to get an estimate of 
the scale heights. To simplify things, representative values of a were determined from the 
average over a range of parameters: D = (0,2), Sgo/Ss,o = (0,2), and Rb/Rq = (0.2,2). 
The average a for disk position R/Rc = 1 was found to be 0.4, meaning the rotation curve is 
on the rising part there, and the average for R/Rc = 4 is 0.04, which is on the fiat part. We 
consider these two cases in what follows. The first case can apply to both the inner regions 
of massive disk galaxies and to most of the visible parts of the disks of dwarf galaxies (e.g.. 



Burstein fc Rubinlll985[ ). 



Figure [TO] shows the profiles of pg{z) and Ps{z) for several dark matter ratios D and 
velocity dispersion ratios S. The solid curves are for a = 0.4 and the dashed curves for 



- 33 - 



a = 0.04; they hardly differ from each other. For all colors and line types, the lower curves 
are for gas and the upper curves for stars. When S = o-g/ag = 1 (blue curves), the gas scale 
height is the same as the stellar scale height and the profiles differ only in amplitude by the 
ratio of column densities (Sgo/Ss,o = 0.2). When S is small, the gas disk is thin and the 
dark matter ratio does not affect the gas scale height (red and green curves for gas are about 
the same). The stellar profile does not depend much on S (blue and green curves for stars 
are about the same), but it does depend on the dark matter ratio in the sense that less dark 
matter makes the stellar disk thicker (red curve for stars is higher than the blue or green 
curves). In all cases, Qs = 2, Vrot = i^max = 8, and Sg_o/Ss,o = 0.2. 

Figure [11] shows the dimensionless scale heights versus the gas fraction in the disk, 
Sgo/Ss,0) on the left, versus the velocity dispersion ratio S in the middle, and versus Qs 
on the right. These are the parameters that cause the scale heights to vary most. In all 
cases, Vrot = "Omax = 8, and a = 0.4. Equations (lAllI) and (\A12\} show why the gas scale 
height increases with S and both scale heights increase with Qs for this normalization. In the 
left and middle panels, Qs = 2. Both dimensionless scale heights decrease with increasing 
Sg o/Es_o because then the gas density increases at fixed Ss^o and this pulls in the star layer. 
Analogous results for a = 0.04 and for a wide range of -Orot are nearly identical to those in 
Figure [TTl 

Representative values for normal spiral galaxy disks in the local universe, where S ~ 
0.2 — 0.4, and Sg o/Ss,o ~ 0.2, are Hg ~ 0.3 and ~ 1-3. For clumpy galaxies at high 
redshift and local dwarfs, S ~ 0.8 — 1, and Sg o/Ss,o ~ so Hg ~ 0.7 and ^s ~ 0.9. The 
ratio of the rotation speed to the velocity dispersion of the stars, the slope of the rotation 
curve, and the ratio of the dark matter to the disk mass do not influence these normalized 
scale heights significantly. 



B. Thickness Correction Factors 

Equation flTI?]) contains an approximation for disk thickness that can be evaluated using 
the full solutions for the vertical density profiles just calculated. If p{z) = Pg{z) + Ps{z) is 
the vertical distribution of total density in a plane-parallel layer, and if this total density 
has a cos{kx) dependence in the direction of the perturbation, as assumed above, then the 
acceleration measured at position xq in the midplane in this direction is 

9(.„) = 2G r d. r d. ^'^>rf"%;°' . (Bl) 

J-oo J^oo Z^+{X- XoY 

This comes from the acceleration 2Gp/r toward a line of mass with mass per unit length 
/i = pdxdz measured at a distance r = {z^ + [x — xq]^)^/^. Converting equation f lBl|) to 
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dimensionless units as above, separating the gaseous and stellar parts of p, setting g = 
—d(j)/dx, and using an equation like f|T9|) with arbitrary correction factors C instead of 
(1 + kHy\ we get 



Note that the cosine term inside the integral gets converted to a sine term after integration, 
and this phase dependence is divided out by the sine in the denominator. This is consistent 
with the fact that a cosine variation in density produces a sine variation in gravitational 
acceleration, and this produces a cosine variation in the potential, which is used in equation 
f lT9|) . The integral was evaluated for a variety of cases over the range q = (0,20) and the 
coefficients Cj = (1 + qHi)^^ were found to be good approximations. 

Figure [12] shows the ratio of the right hand side of equation flB2p to the left hand side 
using Ci = (1 + qHi)^^, plotted as a function of dimensionless wavenumber q. This ratio is 
the correction factor that should be applied to the right-hand side of equation (fT9|) . or to the 
column densities, to make the thickness correction more accurate when the approximation 
(1 + kH)~^ is used. Alternatively, the inverse of this correction factor can be applied to 
the stability parameter Q to get a more accurate stability parameter. Fortunately, the 
correction factor is small, less than 15% in all cases, and only that large near kH ~ 1 — 4. 
In the range kH ~ 1 — 10, the average correction in Figure [12] is 12%, which means that the 
true acceleration that drives the instability is higher than the estimate using the potential of 
equation f[T9|) by an average of ~ 12%. The instability growth rate scales as the square root 
of the acceleration per unit length, so the likely correction to the growth rate is an increase by 
~ 6%. Considering the other approximations used here for interstellar and stellar dynamics, 
this correction to the growth rate is not important. We therefore evaluate the dispersion 
relation using the original Vandervoort approximation, as given by equation ( !T9|) . 




(B2) 
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0246802468 10 
Wavenumber, q Wavenumber, q 

Fig. 1. — The reduction factor T for gravitational force is plotted versus the dimensionless 
wavenumber, g, for five different dimensionless rates, s, indicated by color. Oscillating 
solutions are shown on the left and growing solutions are shown on the right. Dashed lines 
are approximate solutions from equations (6) and (13). 
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Fig. 2. — Dispersion relation for the planar disk instability with only gas and no rotation, 
showing the effect of turbulent dissipation. 6 is the relative energy dissipation rate compared 
to the perturbation crossing rate. Without turbulent dissipation (5 = 0), the dispersion re- 
lation has the usual upper limit for instability at the Jeans wavenumber k = kj. With 
turbulent dissipation, the instability also happens in smaller regions because shocks remove 
energy and promote local gravitational instabilities. Solid lines are for a velocity disper- 
sion that is independent of the size of the perturbation, as in the rest of this paper, and 
dashed lines are for a velocity dispersion that increases as the square root of the size of the 
perturbation. 
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Fig. 3. — Solutions to the dispersion relation between growth rate and wavenumber are 
shown for a pure gas disk on the left and a pure star disk on the right. The different curves 
of the same type and color in the left panel are for different dimensionless dissipation rates. 
High dissipation leads to greater instability at high wavenumber. On the right, the four blue 
curves correspond to four values of Qs including thickness corrections. The dashed green 
curve is for a zero thickness disk. The dotted red curve is a solution for a zero-thickness 
disk that uses a reduction factor J-" evaluated in the case where the dispersion relation is 
oscillatory; such oscillations are unphysical for the assumed low Qs, because these would give 
growth solutions. In comparison to the dashed green curve, the dotted red curve shows the 
importance of using the appropriate expression for the reduction factor. 
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Fig. 4. — Dispersion relations for pure gas disks are shown in cases formerly thought to be 
stable, Qg > 1- The destabilizing effect of the dissipation term involving 6 can be seen. 
Thickness effects are included, and the stability at low q is from the dilution of in-plane 
self-gravity by disk thickness. 
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Fig. 5. — Solutions to the combined gas+star dispersion relation are shown. The left hand 
panel has several two-component cases that are effectively one-component, as a check on the 
equations. The dashed blue curves are two-component with Qs = I, Qg = 0.5, and S = 0.5, 
and with dissipation parameter S = 0.5 on the top and 5 = on the bottom. The middle 
panel shows a range of values for Qs and Qg at fixed S' = 0.5 and 6 = 0.5. The sequence of 
blue curves with decreasing height has Qs = 0.5, 1, 1.5, and 2 for Qg = 1, and the sequence 
of red curves with decreasing height has Qg = 0.5, 1, and 1.5 for Qs = 1. The right hand 
panel shows the effect of dissipation as a sequence of blue curves {6 = 0, 0.1, 0.3, 0.5, and 1 
with 5* = 0.5 from bottom to top) and the effect of the velocity dispersion ratio as a sequence 
of red curves {S = 0.1, 0.2, 0.3, 0.5, and 1 with 6 = 0.5 from right to left). On the right, all 
have Qs = Qg = 1. 
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Fig. 6. — The ratio of perturbed stars to gas (top) and this same ratio normahzed to the 
initial star-to-gas ratio (bottom) are shown versus wavenumber. The order of the panels, 
left to right, and the curve colors and types are the same as in Figure 7. 
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Fig. 7. — Histograms are shown of the normahze ratio of perturbed stars to gas at the 
wavenumber of peak growth. Each histogram is from a series of 800 cases with a range of 
Qs and Qg. Each panel has a different dissipation parameter and velocity dispersion ratio, 
as indicated. As dissipation increases in the high gas dispersion case (from top to bottom 
on the left), the gas plays a more prominent role in the growth of the instability (histograms 
shift to the left). 
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Fig. 8. — The peak growth rate of the gas+star dispersion relation is plotted as a function 
of 1/Qtot for the Wang & Silk (1994) definition, on the left, and the Romeo & Wiegert 
(2011) definition in the middle. In these two cases, the velocity dispersion ratio is 0.2. On 
the right, the velocity dispersion ratio is unity and the two definitions of Qtot are the same 
because the disk is effectively single- component. Each panel contains 20 blue curves and 
20 red curves. Within each curve, Qg varies, while the different curves are for different Qs- 
In the top four rows, blue curves use an expression for Qtot in which thickness effects are 
not considered, even though they are considered in the dispersion relation; red curves use 
expressions for Qtot that include thickness corrections, and are therefore self-consistent. The 
top four panels differ in their value of the dissipation parameter S. The bottom row solves 
the dispersion relation for a disk with zero thickness, which makes the red curves the same 
as the blue curves. The top four rows show how increasing dissipation always lowers 1/Qtot 
at the threshold for instability, Speak = 0. 
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Fig. 9. — These curves plot the average values of l/Qtot where the peak growth rate, Speak, 
is between 0.03 and 03, from Figure 10. These are essentially the stability thresholds for 
the various definitions of Qtot- Blue curves ignore thickness corrections in the definition of 
Qtot, as in Figure 10, while red curves explicitly account for thickness in Qtot (all curves have 
thickness corrections in the actual dispersion relation). The three symbols on the ordinate 
axis are cases with zero-thickness disks. Dissipation lowers 1 / Qtot in all cases. For reasonable 
values of the parameters, Qtot has to exceed 2 or 3 for a dissipative disk to become stable, 
instead of the usual threshold of ~ 1. 
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Fig. 10. — Vertical profiles for gas (quickly falling curves) and stars for two values of the 
dimensionless dark-matter parameter D (defined above equation AlO) and two values of the 
ratio of gaseous to stellar velocity dispersion, S. The abscissa is the dimensionless height 
above the midplane. Solid curves are for a steep rotation curve, a = 0.4, and dashed curves 
are for a nearly flat rotation curve, a = 0.04 (defined below equation A8). Density is 
normalized to k^/ {2tiGQs) (see equations All and A12). These solutions assume = 2, 
Vrot = 8, and Sg,o/Ss,o = 0.2. 
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Fig. 11. — Dimensionless scale heights are plotted versus the gas fraction in the disk on the 
left, versus the gas-to-star velocity dispersion ratio in the middle, and versus Qg on the right. 
Parameters associated with the various colors are indicated at the top of each panel. In the 
left and middle panels, Qs = 2. 
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Fig. 12. — The curves show the correction factors for the in-plane acceleration in a thick 
disk compared to the approximate value using the term (1 + kH). The correction factor for 
the growth rate is approximately equal to the square root of these acceleration factors. The 
correction factors are small, indicating that the usual approximation for disk thickness is 
good. 



